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Abstract 

Downward going atmospheric tau neutrino flux is estimated in the presence of three neutrino 
flavor oscillations for 1 GeV < E < 10 3 GeV. The relative differences between the three and purely 
two neutrino flavor oscillations are elaborated. As an implication, the downward going atmospheric 
tau neutrino flux is compared with the galactic plane tau neutrino flux that is also estimated in the 
presence of three neutrino flavor oscillations. It is pointed out that the galactic plane tau neutrino 
flux dominates over the downward going atmospheric tau neutrino flux until E ~ 10 GeV. 
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I. INTRODUCTION 



In order to perform a meaningful search of extra-atmospheric neutrinos, the atmospheric 
neutrino background is needed to be known under various possible approximations The 
search for astrophysical tau neutrinos is of special interest since a positive search shall not 
only yield the useful information from the surrounding cosmos but shall also corroborate 
the neutrino flavor mixing hypothesis Q|. 

In view of the presently available restricted ranges of the various neutrino mixing param- 



eters in the three neutrino flavor approximation 



as well as the recent relevant detector 



developments 0, it is of some interest to estimate the resulting new atmospheric tau flux. 
Such studies may also provide some examples of new windows to study the cosmos around 
us, that shall open up under the hypothesis of the neutrino oscillations. 

Previously, there is no attempt to estimate the atmospheric tau neutrino flux in the 
context of three neutrino flavor oscillations framework. So far, the atmospheric tau neutrino 
flux is estimated in two neutrino flavor approximation only p|. Recently, it is pointed out 
that the tau neutrino astronomy in the multi-GeV energy range is possible, at least in 
principle. It is so, because the atmospheric v T flux is generally suppressed compared to the 
atmospheric flux, as a result the prospective observation of astronomical v T suffers much 
less background than in the case. The zenith angle (£) dependence of the atmospheric 
background tau neutrino flux is also crucial in determining the possible future prospects 
f|. This is an example of the availability of a new cosmic horizon because of the neutrino 
oscillations, and is true for tau neutrinos only among the three neutrino flavors. 

It is also pointed out that for 0° < £ < 60°, the downward going atmospheric tau 
neutrino flux is minimum, thus providing an opportunity to search for extra-atmospheric 
tau neutrinos in this zenith angle range, in the presence of neutrino oscillations [7]. The 
neutrino oscillation effects are minimal for the down ward going atmospheric neutrinos. As 
a result, the oscillated atmospheric tau neutrino flux provides a minimal background for 
prospective astronomical tau neutrino flux searches. 

Collectively, in these papers it is pointed out that contrary to general expectations, the 
atmospheric neutrino flux does not smear out the astronomical neutrino flux on the Earth 
at low energy (E < 10 3 GeV) once the flavor composition of the incoming neutrino flux is 
taken into account. 
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In this paper, we estimate the atmospheric tau neutrino flux in three neutrino flavor 
mixing approximation. This is an extension of the above mentioned previous works in this 
direction. As an implication of this estimate, we compare it with galactic-plane tau neutrino 
flux in the presence of neutrino oscillations. It is pointed out that the dominance of galactic 
plane tau neutrino flux over the downward atmospheric tau neutrino flux still persists in the 
multi-GeV energy range. 

The considered energy range (1 GeV < E < 10 3 GeV) is interestingly within the reach 



a. 



Tau neutrino 



of presently operating detectors as well as those under active planning 
flavor discrimination is an interesting challenge for the existing/forthcoming astrophysical 
neutrino telescopes. 

This paper is organized as follows. In section II, the electron, muon and tau neutrino 
flux originating from the Earth atmosphere is briefly discussed. In section III, the neutrino 
oscillation effects are studied for these. In section IV, an implication of our estimate for the 
atmospheric tau neutrinos is mentioned, whereas in section V, conclusions are presented. 

II. ATMOSPHERIC NEUTRINO FLUX 

Atmospheric neutrino flux arises when the incoming cosmic rays interact with the air 
nuclei A, in the Earth atmosphere jsj. For 1 GeV < E < 10 3 GeV, the ir, the K productions 
and their direct and indirect decays are the main sources of electron and muon neutrinos, 
both being in the region of conventional neutrino production. The absolute normalization of 
the conventional atmospheric neutrino flux is presently known to be no better than (20—25)% 
in the above energy range |l| . For present estimates, the non tau atmospheric neutrino flux 
is taken from Ref. ^(|. These are atmospheric neutrino flux calculations in one dimension 
without geomagnetic field effects. At higher energy the prompt muon neutrino production 
from .D's dominates over the conventional one [111 ]. 

The atmospheric tau neutrino flux arises mainly from production and the decay of D$ and 
is calculated in Ref. 12, Q]. The Quark Gluon String Model (QGSM) is used in Ref. [k| 
to model the pA interactions. The low energy atmospheric tau neutrino flux is essentially 
isotropic 



12j. For E < 10 2 GeV, the atmospheric tau neutrino flux is obtained by following 
the procedure given in Ref. HQ and re-sca.ing w.,t new cosmic ray flux spectrum, taking 
it to be predominantly the protons [14J. 
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For illustrative purpose, the three downward going (cos £ = 1) atmospheric neutrino fluxes 
F®(E) = d 2 iV°/d(logio-E)d£, in units of cm _2 s _1 sr _1 , estimated using the above description, 
are shown in Fig. Q as a function of the neutrino energy. 

III. NEUTRINO OSCILLATION ANALYSIS: THE ATMOSPHERIC TAU NEU- 
TRINO FLUX 



We shall perform here the three neutrino flavor oscillation analysis. In the context of 



mixing in three flavors, the neutrino mixing parameters are : 5m 2 2 , 5m| 3 , Sm\ 3 , 612, 623, 



13 



and the CP violating phase 5. The 5m 



where i,j = 1,2,3, is the absolute 



difference of the mass squared value of the mass eigenstates and the three 6i/s are the mixing 
angles. 



Presently, the 8m\ 2 an d 0i ? are mainly determined from the data analysis of the solar neu- 
trino flux measurements [la]. The neutrino mixing parameters inferred from solar neutrino 



flux studies are recently confirmed by a terrestrial reactor experiment, namely Kamiokande 



Liquid-scintillation Anti-Neutrino Detector (KamLAND) lldl . The 5m 2 3 and #23 are mainly 
determined from the study of the atmospheric neutrinos 17[. The above neutrino mixing 
parameters inferred from atmospheric neutrino flux studies are also recently confirmed by a 
terrestrial accelerator experiment, namely KEK to Kamiokande long-baseline neutrino ex- 
periment (K2K) The #13 is mainly constrained by the CHOOZ experiment 3- For 
some recent feasibility studies to determine more precisely the #13 and/or 5, see |20j]. The 
total range of 5m 2 irrespective of neutrino flavor in the context of three neutrino flavors is 
1(T 5 eV 2 < 5m 2 < 1(T 3 eV 2 . 

Under the assumption that 5 = 0, the 3x3 Maki Nakagawa Sakata (MNS) mixing matrix 
J in standard parameterization connecting the neutrino mass and flavor eigenstates reads 



2l|: 



/ 



U 



\ 



(1) 



C12C13 S12C13 S13 

-S12C23 - C12S23S13 C12C23 - S12S23S13 S23C13 
S12S23 - C12C23S13 -C12S23 - S12C23S13 C23C13 

where = cos#y and s^- = sin^. The presently available empirical constraints for the 



V 



/ 
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various neutrino mixing parameters give the elements of the above U matrix as 

/ nc/i nK/i n 1 n \ 



U e i U e 2 U e s ^ 



u 



UfjX U^2 U^3 
UtI Ut2 U t s 



\ 



0.84 0.54 0.10 
-0.44 0.56 0.71 
0.32 -0.63 0.71 



(2) 



A recent short summary of the model predictions for the expected magnitude of \U P z\ i s 



given in Ref. 



221 ] . Using Eq. ([TJ, the neutrino flavor oscillation probability formula is \2c 



P(Vn 



U 



E, = P al3 = J2 U 2 ai U% + J2 U m U pi U aj U P3 cos (^) , 

i=l & ^ ij ' 



(3) 



where a, (3 — e, //, r and Ljj ~ AE/Smfj is the neutrino oscillation length. The L in Eq. (J3J) 
is the neutrino flight length. In the Earth atmosphere, it can be estimated using the relation 



L = ^(h 2 + 2R®h) + (i? e cosO 2 - R® cos^, (4) 
where ^ is the zenith angle. The L is essentially the distance between the detector and the 
height at which the atmospheric neutrinos are produced. The ~ 6.4 ■ 10 3 km is the 
Earth radius, and for instance, = 15 km is the mean altitude at which the atmospheric 
muon neutrinos are produced. In general, h is not only a function of the zenith angle £, the 



24] 



neutrino flavor but also the neutrino energy 

We consider the zenith angle range between 0° and 60°. Because in this energy range 
i) as pointed out in Ref. |7J, the atmospheric neutrino oscillation effect is minimal, ii) 
the Earth curvature effects can be ignored, and iii) more relevant for present discussion 
is that in this zenith angle range, the matter effects in neutrino oscillations are negligible 
as the Earth atmosphere matter density traversed by neutrino satisfies: p < p res where 
p res = m N Acos26/y/2G F} with A = Sm 2 /2E ranges between 10~ 17 eV and 10~ 12 eV, for 
the entire considered 5m? and E ranges. Fig. El shows an example of P^ T using Eq. (JBJ) 
for four different £ values where £ = 0° corresponds to the downward going atmospheric 
neutrinos. A relevant remark useful for what follows next is that for E < 20 GeV, the P MT 
is larger by a factor of up to ~ 8 for larger £, owing to the fact that for larger £, neutrinos 
traverse larger distance in Earth atmosphere. 

In the limit #12, #13 — » 0, we obtain the commonly used two neutrino flavor oscillation 
probability formula for Piy^ — ► u T ;E, £) as 

sin 2 2# 23 ■ sin^ I 1.27 



P 



2 / , ^m 2 3 (eV 2 )L(km) 



fir 



E(GeV) 



(5) 



where sin 26*23 = 2C23S23. The neutrino flux F Va (E,£), arriving at the detector, in the 
presence of neutrino oscillations is estimated using 

F Ua (E,0 = Y,P^(E,OF°^(E,0, (6) 

P 

where F® (E, £) are taken according to discussion in section II. The P a p(E,£) is a 3x3 matrix 
obtainable using Eq. (j3J). The unitarity conditions such as 1 — P ee (E,^) = P eft (E,£) + 
P er (E,^) are implemented at each E and £ at which these are evaluated. Fig. El shows 
comparison of the atmospheric tau neutrino flux for £ = 0° and £ = 60° in three and two 
neutrino flavor mixings as a function of the neutrino energy. Note the difference in the 
slope in three and two neutrino flavor approximations because of the presence of second 
Sm 2 scale in three flavor mixing relative to only one Sm 2 scale in purely two neutrino flavor 
approximation. The atmospheric tau neutrino flux for larger zenith angles is larger relative 
to smaller zenith angles for E < 20 GeV (see Fig. 0also). 

The downward going atmospheric tau neutrino flux in the presence of three neutrino 
flavor oscillations can be parameterized for 1 GeV < E < 10 2 GeV as 

F Vt (E)=A- Bx- 1 + Cx- 2 + Dx- 3 - Ex~\ (7) 

where the coefficients are A = 4.20865 ■ 10 -11 , B = 1.57226 ■ 10~ 8 , C = 5.37912 • lO" 6 , 
D = 3.97309 • lO" 5 and E = 2.21708 • 10" 5 with x = E/GeV. The F Vr (E) is in units of 

—2 —1 —1 

cm s sr 



IV. AN IMPLICATION: COMPARISON WITH THE GALACTIC PLANE TAU 
NEUTRINO FLUX 

In this section, we shall estimate the galactic plane tau neutrino flux with three neutrino 
flavor mixing. This is to illustrate a new opportunity we may have to study cosmos having 
estimated the related atmospheric tau neutrino background. 

The galactic plane electron and muon neutrino flux is calculated in Ref. I25J|, whereas 

n 

the galactic plane tau neutrino flux is calculated in Ref. [13J . These calculations consider pp 
interactions inside the galaxy with target proton number density ~ 1/ cm 3 along the galactic 
plane, under the assumption that the cosmic ray flux spectrum in the galaxy is constant at 
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its locally observed value. We emphasize that because of above uncertainties, the galactic 
plane neutrino fluxes should only be considered as reference upper limits. 

Following Ref. Q], the galactic plane neutrino flux for 1 GeV < E < 10 3 GeV is obtained 
by re-scaling w.r.t new cosmic ray flux spectrum. The tau neutrino production is rather 
suppressed in the galactic plane relative to muon neutrino production. In general, the non 
tau neutrino flux is larger than the tau neutrino flux for E < 10 3 GeV from the above two 
sites. 

In order to estimate the oscillation effect, the distance L for galactic plane neutrinos can 
be taken as ~ 5 kpc, where 1 pc ~ 3 x 10 13 km. This imply that Lij <C L, namely the 
galactic plane neutrino flux oscillate before reaching the earth. Thus, the galactic plane 
neutrino flux is averaged out for the whole range of 5m 2 in the entire considered energy 
range. Under the assumption of averaging over rapid oscillations, we obtain from Eq. (jSJ) 

3 

i=i 



The 3 x 3 P matrix using Eq. (j2J) in above Eq. is 



/ 



P = 



P P P 

J- p.p. J- Pti. -*- P: 



\ 



P P P 

\ P P P I 

\ 1 er 1 fir A tt J 



( 0.59 0.23 0.20^ 



0.23 0.39 0.39 
0.20 0.39 0.42 



(9) 



Note that this P matrix is independent of not only 5m 2 but also E (an d £). A rather similar 
3 x 3 P matrix was obtained for both vanishing 5 and #13 in Ref. |26(. 

Using Eq. © and Eq. Q, the galactic plane tau neutrino flux is estimated in the presence 
of three neutrino mixing. It is then compared with the atmospheric tau neutrino flux in Fig. 
0] The comparison includes the two neutrino flavor approximation effect with maximal 23 
mixing also. From the figure, it can be seen that the galactic plane/non- atmospheric tau 
neutrino flux starts dominating over the downward going atmospheric tau neutrino flux even 
for E as low as 10 GeV in the presence of neutrino oscillations, depending upon the incident 
zenith angle for Atmospheric neutrinos (the galactic plane neutrino flux is independent of 
the zenith angle). This is a very specific feature of tau neutrinos, and is absent for electron 
and muon neutrinos. This specific behavior has to do with the neutrino oscillations. 
The galactic plane tau neutrino flux in three neutrino mixing is smaller by a factor ~ 2/3 
relative to two neutrino flavor mixing, irrespective of the neutrino energy. 



Fig. E] indicates that zenith angle dependence of the total tau neutrino flux can at least 
in principle help to distinguish between atmospheric and non-atmospheric tau neutrino flux. 
The galactic plane tau neutrino flux dominance is independent of number of oscillating 
neutrino flavors (two or three). The galactic tau neutrino flux transverse to the galactic 

n 

plane is three to four orders of magnitude smaller than the galactic plane one |13l |. 

Fig. |3] shows the three downward going atmospheric neutrino fluxes in the presence of 
three neutrino flavor mixing, along with the corresponding galactic plane neutrino fluxes. 



Here the contribution from the charm production is also taken into account |27[. The 
galactic plane three neutrino flavor fluxes are approximately equal because of Eq. (jHJ). The 
unique behavior of atmospheric oscillated tau neutrino flux relative to electron and muon 
neutrino flux is evident. The oscillated galactic plane tau neutrino flux thus dominates 
until approximately 10 GeV over the corresponding atmospheric neutrino background. On 
the other hand, the same occurs at E ~ 10 5 GeV and E ~ 10 6 GeV for electron and muon 
neutrino fluxes in the presence of three neutrino flavor mixing from the two astrophysical 
sources under discussion. The change in the slope of the atmospheric electron and muon 
neutrino fluxes is a reflection of the corresponding change in the slope of cosmic-ray flux. 

In general, the three downward atmospheric neutrino fluxes provide an energy 
dependent background to the incoming astronomical neutrino flux in the presence of neu- 
trino oscillations. This observation may have some relevance for the forthcoming neutrino 
telescopes with the prospective neutrino flavor discrimination capabilities. 

The galactic plane tau neutrino flux for 1 GeV < E < 10 3 GeV in the presence of three 
neutrino flavor mixing can be parameterized as 

r -I -2.74 

F Ut (E) = 4.38 ■ 10~ 6 ■ E im E + 2.15 exp(— 0.21 v^E) , (10) 

where F Ut (E) is in units of cm~ 2 s -1 sr -1 and on r.h.s. E is in units of GeV. 

We have estimated the galactic plane tau neutrino induced shower event rate for the 

I I 

forthcoming Mega ton (Mt) class of detectors |28j, to indicate the limited prospects offered 
by below TeV (1 TeV = 10 3 GeV) tau neutrino astronomy to search for extra- atmospheric 
astrophysical neutrino sources in this energy range, as envisaged presently. The tau neutrino 
event rate iV even t for 10 GeV < E < 10 3 GeV can be approximately estimated as 

IVevent = / F Vt (E) ■ OV T — S hower {E) — . (11) 
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The F Vt (E) is the galactic plane tau neutrino flux in the presence of neutrino oscillations in 
units of cm _2 s _1 sr _1 and is given by Eq. (JTUJ). The <j Ut shower (E) is a^° N (E)/Mt, where one 
megaton of water contains ~ 5.5 • 10 35 nucleons j^. The (T UT ^^n Wer (E) is taken from the 
web site of Cern Neutrinos to Gran Sasso (CNGS) experiment |30j]. The N event is found to 
be in the range 1 to 10 in units of Mt _1 yr _1 over 2n ■ sr with a 3 to 5 years of data taking. 



V. CONCLUSIONS 

Atmospheric tau neutrino flux is estimated for 1 GeV < E < 10 3 GeV in the presence of 
neutrino oscillations. The three neutrino flavor oscillation study is carried out. The relative 
difference between two and three neutrino flavor oscillation analysis is elaborated. 

As an implication of this study, the atmospheric tau neutrino flux is compared with 
the galactic plane tau neutrino flux with three neutrino mixing. It is pointed out that 
the dominance of galactic plane tau neutrino flux persists over the downward going (£ = 
0°) atmospheric tau neutrino flux for E > 10 GeV relative to E > 8 GeV in two flavor 
approximation. The same dominance shifts to E > 10 GeV relative to E > 9 GeV for 
f = 60°. 

In the multi-GeV energy range, the dominance of the galactic plane tau neutrino flux 
over the atmospheric tau neutrino flux for incident zenith angle between 0° and 60° in the 
presence of neutrino oscillation both in two and three neutrino mixing is unique to tau 
neutrinos only. The galactic electron and muon neutrino flux dominance starts only for 
E > 10 5 GeV irrespective of mixed (two or three) neutrino flavors. 

Summarizing, a possibility of a neutrino flavor dependent new window in low energy 
neutrino astronomy seems to exist. Our present study is an attempt to justify this existence. 
Namely, for the incident neutrino energy range between 10s of GeV up to 1000s of GeV, 
in which the presently data taking large neutrino detectors such as the Antarctic Muon 



And Neutrino Detector Array (AMANDA) experiment 



experiment 



at the south pole and the Baikal 



32| in the lake Baikal in Russia are not (yet) optimized to perform extra- 
atmospheric neutrino astronomy, the upcoming Mega ton class of detectors (or even the 
currently existing ones) should at least in principle be able to do possibly a meaningful 
neutrino astronomy based on prospective neutrino flavor identification. 
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FIG. 1: The down ward going atmospheric electron, muon and tau neutrino fluxes as a function 
of the neutrino energy. Absence of neutrino oscillations is assumed here. More details are given in 
the text. 
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P^ x for four different zenith angles 
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FIG. 2: An illustration of P^ T for four different £ values, as a function of the neutrino energy. 
Three neutrino flavor mixing is assumed here. 
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Atmospheric Tau Neutrino Flux 
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FIG. 3: The downward going atmospheric tau neutrino flux in the three and the two neutrino 
mixings for f = and f = 60 as a function of the neutrino energy. 
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FIG. 4: Left panel: The comparison of the galactic plane and the downward going atmospheric 
tau neutrino flux for £ = 0° as a function of the neutrino energy. Three flavor neutrino mixing is 
assumed. Also shown is two neutrino flavor mixing results with maximal 23 mixing. Right panel: 
Same as left panel but for £ = 60°. 
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Atmospheric and Galactic Neutrino Fluxes 
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FIG. 5: The down ward going atmospheric electron, muon and tau neutrino fluxes as a function 
of the neutrino energy in the presence of three neutrino flavor mixing. The galactic plane neutrino 
flux estimated in three neutrino neutrino flavor mixing is also shown for comparison. 
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